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We study robustness of genuine multipartite entanglement for a system of three qubits under 
collective dephasing. Using a computable entanglement monotone for multipartite systems, we 
found that almost every state is quite robust under this type of decoherence. We analyze random 
states and weighted graph states at infinity and found all of them to be genuinely entangled. 
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Quantum entanglement is essential for several practi¬ 
cal applications in the emerging field of quantum Infor¬ 
mation. This feature inspire us to develop a theory of 
entanglement and several methods to partially avoid the 
unalterable process of decoherence [U, Q • As we desire to 
generate and manipulate entanglement in experiments, 
therefore it is essential to study the effects of various envi¬ 
ronments on entanglement Q . This problem has received 
considerable attention for bipartite and multipartite sys¬ 
tems Few studies considered the life time of en¬ 

tanglement under decoherence ill , however life time 
may not characterize the decoherence process appropri¬ 
ately. The life time of entanglement may be long, but 
the actual amount of entanglement might be very small 
after a short time, making the decohered state practi¬ 
cally useless for information processing tasks Q. Some 
authors have studied only l^artite aspects of entangle¬ 
ment of multipartite states Q , nevertheless these studies 
can capture only partial information because genuine en¬ 
tanglement is distinct in nature than the entanglement 
among various partitions. Unless there is an adequate 
theory of multipartite entanglement, one can talk about 
lower bounds of entanglement instead of its exact value 
Q. The precise expression of multipartite entanglement 
was only computed for particular type of decoherence and 
definite quantum states [HI- In addition, in order to bet¬ 
ter judge the dynamical behaviour of a state, we need to 
compare its dynamics with dynamics of random states. 
In our previous studies, we focussed on the dynamics of 
multipartite states by investigating a computable mea¬ 
sure of genuine multipartite entanglement for several de¬ 
coherence models (l7l| . 

In this Letter we extend our studies to examine the 
robustness of three-qubit states as well as random states 
under collective dephasing. Recent progress in the the¬ 
ory of multipartite entanglement has enabled us to study 
decoherence effects on actual multipartite entanglement 
and not on entanglement among bipartitions. In partic¬ 
ular, the ability to compute genuine negativity for mul¬ 
tipartite systems has eased this task [l^ [l^. We find 
that due to collective damping effects, decoherence-free- 
subspace exist in this system. This space is spanned by 
the W-type states. Hence the W-type states are most 
robust being uneffected by decoherence process [l^. Al¬ 
though the entanglement of GHZ states is degraded by 


decoherence process, nevertheless, they are also quite ro¬ 
bust and loose their genuine entanglement only at in¬ 
finity. Interestingly, almost every pure random state is 
robust under collective dephasing. 

We consider the three qubits to be spins for exam¬ 
ple and impose a stochastic magnetic field B{t) acting 
on three qubits together. In the interaction picture, the 
Hamiltonian of the qubits and classical noisy field can be 
written as (with h = 1) 

H{t) = +af+ a'^), (1) 

where fi is the gyromagnetic ratio and are 

the standard Pauli matrices in the computational basis 
{|0), |1)}. The stochastic magnetic field refer to statisti¬ 
cally independent classical Markov process satisfying the 
conditions 

{B{t)Bit')) = ^5{t-t'), 

(H(t))=0, (2) 

where (■ • •) is ensemble time average and T denote the 
phase damping rate for collective dephasing. 

We choose the standard computational basis { |000), 
|001),..., |111)}. The time-dependent density matrix 
for three-qubits can be obtained by taking ensemble av¬ 
erage over the noise field, i.e., p{t) = (pst(Q), where 
Pst(t) = U{t)p{0)W{t) and U{t) = exp[—i dt'H(t')\. 
We assume that there are no initial correlations be¬ 
tween the principal system and stochastic field, that is, 
p(0) = PS ® Pr- There are various ways to compute the 
time evolved density matrix for three qubits, however we 
prefer the master equation approach to solve the system. 

We consider a system ”S” interacting with a reservoir 
”R”. The combined density operator is given as psR- 
The reduced density matrix ps for system alone can be 
obtained by taking trace over reservoir degrees of free¬ 
dom, that is, PS = Tri{(ps/{). In the interaction picture, 
the equation of motion is given as (with /i = I) 

ipsR = [H, psR(,t)]. (3) 

As the interaction between system and reservoir is weak, 
therefore we look solution of the form pssit) = Ps{t) ® 
PR{ti) + Pcit), where for consistency we demand that 
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Tr/{(/9c(t)) = 0. Substituting this form and Hamiltonian 
o, retaining terms upto order , and after doing some 
simplification, we get 

f 1(0-^ + + + '^f + o-f)pW 

- 2 (cr^ + erf + erf) p{t) (erf + erf + erf) 

+ P{t){<^z + 0-f + o-f + o-f + erf) |. (4) 

We have used the relations @ to obtain the above equa¬ 
tion. It is now straight forward to obtain the most general 
solution. Due to collective dephasing, there appear the 
decoherence free subspaces (DFS) for W type states [^ . 
This DFS is spanned by basis { |001), |010), |100), |011), 
1101), 1110)}. Hence the W type states are most robust 
states under this type of decoherence model. 

To present the concept of genuine entanglement, we 
take a system of three parties A, B, and C. A state is 
defined as biseparable if it can be written as a mixture 
of states which are separable with respect to different 
bipartitions, that is 



A state is genuinely multipartite entangled if it is not 
biseparable. This definition can readily be extended to 
A^-partite systems. 

Genuine entanglement can be detected and charac¬ 
terized via the technique of p ositive partial transpose 
mixtures (PPT mixtures) A bipartite state p = 

® l^)(^l is PPT if its partially trans¬ 
posed matrix PjiM |*)(j| ® |^)(^| is positive 

semidefinite. As separable states are always PPT [2l| and 
therefore the set of separable states with respect to some 
partition is contained in a larger set of states which has a 
positive partial transpose for that bipartition. The states 
which are PPT with respect to fixed bipartition may be 
called Pb\ac' Pc\ab- whether a state 

can be written as a mixing of PPT states 

pPPTm.. ^ pPPT^ ^ pPPT^ ^ pPPT^ ^ (g) 

which is called a PPT mixture. Because any biseparable 
state is a PPT mixture, any state which is not a PPT 
mixture must be genuinely entangled. The main ben¬ 
efit of employing PPT mixtures instead of biseparable 
states is that PPT mixtures can be f ully characterized 
by semidefinite programming (SDP) (22|. In addition, 
the set of PPT mixtures serves a very good approxima¬ 
tion to the set of biseparable states and delivers the best 
known separability criteria for many cases, nevertheless 
there are multipartite entangled states which are PPT 
mixtures [l^. A state is a PPT mixture if and only if 

minTr(>Vp) (7) 

has a positive solution under the constraint that for all 
bipartitions M\M 

W = Pm + , with 0 < Pm < 7,0 < Qm < I- (8) 


The constraints reflect that W is a decomposable entan¬ 
glement witness for any bipartition. If this minimum is 
negative then p is not a PPT mixture and hence genuinely 
entangled [l^. Since this is a semidefinite program, the 
minimum can be efhciently computed and the optimality 
of the solution can be certified [2^. We use the pro¬ 
grams YALMIP and SDPT3 [2^ to solve SDP. We also 
use implementation which is freely available [^ . For 
more explanations, see Ref. 0. 

The absolute value of the minimization was shown to 
be an entanglement monotone 18|. This monotone has 
been called genuine negativity 19|. In this work we use 
this measure by denoting it A(pb For bipartite systems, 
it is equivalent to negativity [25j | and for a system of 
qubits, it is bounded by E{p) < 1/2 [^ . 

Two inequivalent genuinely entangled states for three 
qubits are the GHZ states \GHZi) = l/-\/2(|000)-|-|lll)) 
and the W states \W) = 1/V3(|001) -b |010) -b |100)). 
Entanglement monotone for GHZ state has a value of 
E{\GP[Z){GP[Z\) = 1/2, whereas for the W state, its 
value is E{\Wz)(Wz\) ~ 0.443. As the W state lives in 
DFS, therefore it is already robust, hence we only focus 
on GHZ states. Another locally equivalent GHZ state is 
\GHZ 2 ) = l/-\/2(|001) -b |110)). In order to make mean¬ 
ingful comments on robustness of these states, we com¬ 
pare the dynamics with that of random pure states and 
weighted graph states. The method to generate these 
random states is described in detail in Ref. [la. 

We argued before that the life-time of entanglement 
may not lead to conclusive results. In addition, start¬ 
ing with initial states having different amount of entan¬ 
glement, the lifetime of entanglement may not tell the 
relative decay. Moreover, the comparison of the values 
of E[pi{t)] for different t and initial states i is not use¬ 
ful. To overcome these problems, we may consider the 
logarithmic derivative 


pit) 


diHm]) 

dt 


d/dt [E{t)] 

m 


(9) 


where E(t) is the entanglement monotone [HI- This 
quantity describes the relative decay of entanglement, 
and it is now possible to compare states with different 
initial entanglement. We note that for a state where 
the entanglement just decays exponentially, r](t) is con¬ 
stant and the inverse of the half-life. We stress that as 
in our previous study [13 , we use logarithmic derivative 
of genuine entanglement to claim the robustness of all 
states. The term ’’robustness of entanglement” was used 
in Refs. [2ji3> but in it was used as a kind of quan¬ 
tification of entanglement. In [ 2 ^, the authors have iden¬ 
tified most robust states under local decoherence used 
the definition introduced before whereas, the authors 
of Ref. [ 2 ^ considered asymptotic long-time dynamics of 
initial states and identified two different classes of states, 
one class is fragile even there remains some coherence 
in the system and the second class as most robust states 
which become disentangled only when decoherence is per¬ 
fect. However, our approach is different than these stud- 
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FIG. 1. Logarithmic derivative is plotted against parameter 
Ft GHZ states, mean values for weighted graph states and 
random pure states. 

ies. 

Figure [T] shows logarithmic derivative plotted against 
parameter Ft. The constant dotted line at —0.5 is for 
\GHZ 2 ) state whereas constant solid line at —4.5 is for 
\GHZi) state. This means that both states remain gen¬ 
uinely entangled till infinity and decay only exponen¬ 
tially. This claim can also be verified from the fact that 
as all zero elements in the density matrix remain zero, 
therefore GHZ states is genuinely entangled if \p 27 \ = 
g-ri/2/2 > 0 for \GHZ 2 ) state and Ipisl = > 0 

for \GHZi) state [S^l, which is obviously the case here. 
The thick dashed line is the mean value of logarithmic 
derivative for 100 random pure states whereas the thick 
dashed-dotted line is the corresponding mean value for 
100 weighted graph states. We observe that these two 
mean values converge to the value of \GHZ 2 ) state. This 
means that almost every random state is quite robust. 
This interesting feature is due to presence of decoherence- 
free-subspaces (DFS). 

In Figured! we plot logarithmic derivative r]{t) for 100 
random pure states. We observe that every random state 
is robust under collective damping. We also plot the 
mean value of r](t) (thick solid line). From this data we 
also obtain an error estimate to indicate the reliability of 
measure. This can, for instance, be defined as a confi¬ 
dence interval 

GI = p±VS, (10) 

where /i stands for mean value and S for variance of quan¬ 
tity being measured. Note, however, that this is not a 
confidence interval in the mathematical sense. The thick 
dashed lines are the confidence intervals for the mean 
value. 

Figure |3| shows the logarithmic derivative for 100 
weighted graph states plotted against Ft. As in the case 
of random pure states, every state here tends towards 
the value of —0.5. The thick solid line is the mean value 
of r]{t) whereas dashed lines are the corresponding confi¬ 
dence intervals. 

Figure m is the bar plot of entanglement monotone for 
n = 100 random pure states at t = 00 . It can be seen 
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FIG. 2. Logarithmic derivative for random pure states plotted 
against parameter Ft. See text for details. 
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FIG. 3. Logarithmic derivative for weighted graph states plot¬ 
ted against parameter Ft. See text for details. 
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FIG. 4. Entanglement monotone for asymptotic random 
states plotted against parameter n. 


that each state has strictly positive value of entangle¬ 
ment monotone. Hence all random pure states remain 
genuinely entangled even at infinity. 

Figure [5] depicts the bar plot of entanglement mono¬ 
tone for n = 100 weighted graph states at t = 00 . Once 
again all states remain genuinely entangled forever, how¬ 
ever, the amount of entanglement monotone is lower as 
compared with random pure states. 

We analyze the entanglement properties of asymptotic 
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FIG. 5. Entanglement monotone for asymptotic weighted 
graph states plotted against parameter n. 


random states where the complete decoherence has oc- 
cured and 7 = = 0. The non-vanishing density 


matrix elements at t = 00 are {pn, P 22 , P 23 , P25, P32, 
P33, P35, P44, P46, P47, P52, P53, P55, P64, P66, P67, P74, 
P76, P77, Pss}- 

We have studied the robustness of genuine multipar¬ 
tite entanglement under collective dephasing. Due to ap- 
pearence of decoherence free subspaces, our findings show 
that almost every quantum state in this type of dynamics 
is robust. As the W type states live in this DFS, there¬ 
fore they are effectly immune to decoherence and hence 
most robust states. Interestingly, We have compared the 
dynamics of specific states with dynamics of random pure 
states and weighted graph states. Again, we found that 
every single state thus generated is robust and remain 
genuinely entangled throughout the dynamical process. 
This situation is in contrast with two qubits case where 
there are always some fragile quantum states under col¬ 
lective dephasing. We believe that our findings would be 
motivational for an experimental verification. 

The author is grateful to referees for their positive com¬ 
ments and suggestions. 
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